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Abstract 

We present a brane- world scenario in which two regions of AdS§ space-time are glued together 
along a 3-brane with constant positive curvature such that all spatial dimensions form a compact 
manifold of topology S 4 . It turns out that the induced geometry on the brane is given by 
Einstein's static universe. It is possible to achieve an anisotropy of the manifold which allows 
for a huge hierarchy between the size of the extra dimension R and the size of the observable 
universe Rjj at present. This anisotropy is also at the origin of a very peculiar property of our 
model: the physical distance between any two points on the brane is of the order of the size of 
the extra dimension R regardless of their distance measured with the use of the induced metric 
on the brane. In an intermediate distance regime R <C r <C Ru gravity on the brane is shown to 
be effectively 4-dimensional, with corresponding large distance corrections, in complete analogy 
with the Randall-Sundrum II model. For very large distances r ~ Ru we recover gravity in 
Einstein's static universe. However, in contrast to the Randall-Sundrum II model the difference 
in topology has the advantage of giving rise to a geodesically complete space. 
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1 Introduction 



Recent suggestions that large CP-jSl or infinite extra dimensions are not necessarily in conflict 

with present observations provide new opportunities for addressing several outstanding problems 
of modern theoretical physics like the hierarchy problem [J— H3,[S|, |Zj or the cosmological constant 
problem [5]-|lUj. In the course of this development and inspired by string theoretical arguments 
the notion of brane world scenarios emerged in which the usual Standard Model fields are supposed 
to be confined to a so-called 3-brane, a 4-dimensional sub-manifold of some higher-dimensional 
space-time. As shown in [f)] also gravity can appear to be effectively 4-dimensional for a brane- 
bound observer provided the conventional scheme of Kaluza-Klein compactification |12j-|15j is 
replaced by a compactification using non-factorizable (also called warped) geometries (see also (Hj). 
This triggered an immense research activity in theories involving 3-branes with interests ranging 
from elucidating the global space-time structure of brane world scenarios, properties of gravity, 
cosmology and brane cosmological perturbations, generalizations to higher dimensions etc. While 
the possibilities are rich, realistic scenarios remain rare. For example the simple Randall-Sundrum- 
II model |S] faces the problem of being geodesically incomplete [1611171 ITS] . It is therefore reasonable 
to look for adequate alternatives or generalizations to the Randall-Sundrum-II model which avoid 
the above mentioned problems while sharing its pleasant feature of the effective 4-dimensional low 
energy gravity on the brane. 

In this paper we present a 5-dimensional brane-world model which solves the geodesic in- 
completeness of the Randall-Sundrum II model while preserving its phenomenological properties 
concerning the localization of gravity. 

We try to illustrate our motivation for considering a particular geometry by using the simple 
picture of a domain structure in extra dimensions resulting from a spontaneously broken discrete 
symmetry. The associated Higgs-field takes different values in regions separated by a domain wall, 
which restricts the possibilities of combining domain walls depending on the global topology of 
the space-time under consideration. We concentrate on the 5-dimensional case, see [7].|19|-[2"1] for 
higher dimensional constructions. 

Let us analyze a few simple cases: 

• Non-compact extra dimension: y G (— oo, +00). In this case, by choosing the origin y = to 
coincide with the position of the brane, we obtain two non-overlapping regions (— 00, 0) and 
(0, +00) and we can thus consistently have one brane in such a theory. An example for this 
configuration is provided by the Randall-Sundrum II model [S]. 

• Compact extra dimension: y G [0,27r]. We now consider two possibilities, depending on the 
spatial topology of our manifold: 

— (a) M 3 x S : If the ordinary dimensions are supposed to be non-compact, it is not 
possible to consistently put only one brane in the extra dimensions. At least two branes 
are needed. 

— (b) S 4 : If all spatial coordinates are part of a compact manifold, the simple picture 
shown in Fig. ^ seems to suggest that it is possible to consistently put a single brane in 
the bulk space-time. 

To the best of our knowledge case (b) has not yet been considered in the literature and our aim is 
to present such a construction. 
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This paper is organized as follows: in Section [21 we discuss the basic geometric and topological 
properties of our brane- world scenario like Einstein equations, junction conditions and the distance 
hierarchy between the extra dimension R and the observable universe Ry. Section El is dedicated to 
the study of geodesies and the demonstration that our model does not suffer from being geodesically 
incomplete. In section 0] we present a detailed computation of the propagator of a massless scalar 
field in the given background serving as an easy, phenomenological approach to the study of gravity. 
We discuss the behavior of the two-point function in three different distance regimes. It turns 
out that the computations are rather technical and we therefore collect large parts of it in four 
appendices. We eventually draw conclusions in section 03 

Upper Hemisphere 




Figure 1: Compact topology. S^-brane embedded in S A . +1 and —1 symbolically represent the 
different vacuum expectation values taken by the Higgs field in different regions of space. 



2 The background equations 

The aim of this section is to present the topology and geometry of the brane- world model which we 
set out to study but also to motivate how this particular model emerged through imposing several 
physical conditions on the more general set of solutions. 

2.1 Einstein equations 

We would like to construct a space-time with all spatial dimensions being part of the same compact 
manifold. One more motivation for this lies in a possible solution of the strong CP problem within 
theories with extra dimensions |2S]. The overall topology would then be given by R x K, where 
R represents the time-coordinate and K any compact manifold. As announced, we will restrict 
ourselves to the case of one extra dimension and an induced metric on the brane characterized by a 
spatial component of geometry S 3 . The idea is to combine two 5-dimensional regions of space-time 
dominated by a cosmological constant A in such a way that the border of the two regions can be 
identified with a 3-brane, constituting our observable universe. As pointed out in J2E| the manifold 
K has to be highly anisotropic in order to single out the small extra dimension from the three usual 
ones. It is not a priory clear whether the Einstein equations allow for such solutions at all and if 
so whether gravity can be localized on the brane in such a setup. 

We choose the following ansatz for the 5-dimensional metric consistent with the above require- 
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ments: 



ds 2 = gMN<ix M dx N 



-a 2 '(9) dt 2 + Rfj 7 2 (#) dnj + R 2 d6 2 



(2.1) 



where 9 £ [— §, f] denotes the extra coordinate and R and Ru are constants representing the size 
of the extra dimension and the size of the observable universe at present, respectively, dil^ denotes 
the line element of a 3-sphere: 



dm 



Q *1 Q Q O O 

d(p l + sin </?i (i<^ 2 + sin Lp\ sin ^2 d<p 3 



(2.2) 



where ^2 belong to the interval [0, it], 923 to [0, 2ir] . Capital latin letters M, N, .. will range from 
to 4. In order to obtain a compact space, we require 7(±§) = 0. 

Metrics which can locally be put into the form (|2.1|) are sometimes referred to as asymmetrically 
warped metrics due to two different functions (a and 7 in our notations) multiplying the temporal 
and spatial part of the 4-dimensional coordinate differentials. Their relevance in connection with 4- 
dimensional Lorentz- violation at high energies was first pointed out in |27] . More recent discussions 
of this subject can be found in [2B] - pIIIl ° r EH an d references therein. A common prediction of 
theories of this kind is that dispersion relations get modified at high energies. For a field-theoretical 
discussion of Lorentz violating effects in the context of the Standard Model of particle physics 
see |3U]. 

The Einstein equations in 5 dimensions with a bulk cosmological constant A and a stress-energy 
tensor Tmn take the following form: 



Rmn — 2 R 9m n + A qmn = -Tp Tmn 



(2.3) 



where M is the fundamental scale of gravity. We choose to parametrize the stress-energy tensor in 
the following way, consistent with the symmetries of the metric: 



e (0)<0. 



8)e{9), T% = e e (9) 



(2.4) 



where as indicated the above diagonal components depend only on the extra coordinate 9. Lower 
case latin indices i,j label the coordinates on S 3 . 

Using the metric ansatz (|2.1|) together with the stress-energy tensor (|2.4|) the Einstein equations 
()2.3|) become 
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where ' denotes differentiation with respect to 6. The conservation of stress-energy, or equivalently, 
the Bianchi-identities lead to the following constraint relating the three independent components 
eo, e and eg: 



a 7 
- +3^ 

a 7 



a 

e 

a 



3^-e 

7 



0. 



(2. 



x Our conventions for the metric and the Riemann tensor are those of reference I2(il . 
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2.2 Vacuum solution 



In this paper we do not intend to provide a field theoretical model which could generate the geometry 
we are about to describe. Our aim is to study a singular brane, located at Of, = separating two 
vacuum regions in the bulk. Of course the solutions to the Einstein equations in vacuum with a 
cosmological constant term are nothing but the familiar de-Sitter (A > 0), Minkowski (A = 0) and 
anti-de-Sitter (A < 0) space-times. To see how these geometries can be recovered using the line 
element (|2.1j) we set the right hand sides of (|2.5j) - (|2.7j) equal to zero and subtract eq. (|2.7|) from 
eq. 1)2. 5 Jl to obtain a = 07', where c is a constant. Putting this back in eqs. (|2.5j) - (|2,7|) . we are left 
with only one independent differential equation 



_3_ 

R 2 



7/ 7 \RuJ T 



+ A = 0. (2.9) 



The solution of (|2.9|) is trivial and by means of a = cy' we find 

, . sinh \uj (f - 101)1 , cosh \oj (5 - |0|)1 

7W = — L . ,y ' 1)1 , <t(o) = — L ,7 ' 1)1 , 2.10 

smh(u;-|) cosh(w^) 

provided that Rjj = RsiTih(u^)/oj and uj 2 = —AR 2 /6. Notice that the solution (j2.10j) does not 
contain any integration constant because we already imposed the boundary conditions 7(rb?) = 
and 7(0) = <r(0) = 1. Moreover, we chose the whole setup to be symmetric under the transformation 
— > — 0. From (|2.1U|) it is now obvious that locally (for > and 8 < 0) the line element (j2.1|) 
correctly describes de-Sitter, Minkowski and anti-de-Sitter space-times for imaginary, zero and real 
values of uj, respectively. For our purposes only the AdS solution will be of any interest as we will 
see in the following section. A simple change of coordinates starting from 1)2. 1|) and (|2,10|) shows 
that the AdS-r&dius in our notations is given by RacIS = R/u = \/~ 6/A. 2 At this stage, the 
validity of the vacuum solution (|2.1Uf) is restricted to the bulk, since it is not even differentiable 
in the classical sense at 8 = 0. In order to give sense to (|2.1()|) for all values of 8 we will have to 
allow for some singular distribution of stress-energy at 8 = and solve (|2.5f) - (|2.7|) in the sense of 
distributions. 



2.3 The complete solution for a singular brane 

As announced, we now refine our ansatz for the energy momentum tensor ()2.4j) to allow for a 
solution of eqs. (|2,5j) - (|2.7|) in the whole interval —n/2 < 8 < tt/2: 

eoW ^ M, <9) _ e m, M-cm. (2 .n) 

After replacing the above components ()2.11j) in eqs. ()2.5|) - (|2,7|) and integrating over 8 from —r/ to 
77, followed by the limit jj->0we find: 

co, (2.12) 
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2 See also footnote 
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where the symbol 



is used to denote the jump of a quantity across the brane defined by: 



:/: ElimJ/^-JH)]. (2.15) 
77 — >o 



Note that in the above step we made use of the identity 

(2.16) 

7 V7/ \7/ 

together with the continuity of 7 on the brane. Specifying (|2.12j) and Q2.13[) to (|2.10|) we have: 

co = -|-M 3 | coth , (2.17) 



c =-— M 3 - 

47T R 



tanh ( u)—) + 2 coth f co— 
2 J V 2 



(2.18) 

c e = 0. (2.19) 

Eqs. ()2.17j) and (|2.18|) relate the energy-density and the pressure of the singular brane to the bulk 
cosmological constant A (via u>), the size of the extra dimension R and the fundamental scale of 
gravity M. With the above relations (|2,17|) - (|2.19|) we can now interpret 1)2.10(1 as a solution to 
the Einstein equations (|2.5jl - (|2.7|) in the sense of distributions. Also the stress-energy conservation 
constraint (|2.8J) is satisfied based on the identity 5(x) sign(x) = again in the distributional sense. 

Note that cq/c — > 1 in the limit u — ► 00. Moreover, for larger and larger values of u, cq and c 
approach the brane tension of the Randall- Sundrum II model and the above eqs. (|2.17|) and (|2.18|) 
merge to the equivalent relation in the Randall- Sundrum II case. This is no surprise since taking 
the limit u> — > 00 corresponds to inflating and flattening the 3-brane so that we expect to recover 
the case of the flat Randall-Sundrum II brane. 

We finish this section by the discussion of some physical properties of our manifold. We first 
observe that its spatial part is homeomorphic to a 4-sphere S 4 . This is obvious from the metric 
(|2.1|) and the explicit expression for 7 given in (|2.10j) . Geometrically, however, our manifold differs 
from S 4 due to the high anisotropy related to the smallness of the extra dimension. The ratio of 
typical distance scales in the bulk and on the brane is given by 

R u) , 

= 1 — r- 2.20 

Ru sinh(^f) V 7 

It is now clear that the above ratio (|2.20|) can only be made very small in the case of real uj (AdS- 
space-time). For the size of the observable universe we take the lower bound Ru > 4 Gpc ~ 10 28 cm 
while the size of the extra dimension is limited from above j^]: R < 10~ 2 cm, leaving us with 
u > 50. 

Finally we would like to point out a very peculiar property of the manifold under consideration: 
as it can immediately be deduced from the line element (|2.1j) and (|2,10|) . any two points on the 
brane are separated by not more than a distance of the order of R regardless of their distance as 
measured by an observer on the brane using the induced metric. 



3 Geodesies 

It is well known that the Randall Sundrum-II model is timelike and lightlike geodesically incom- 
plete ^1 El CHI which means that there exists inextendible timelike and lightlike geodesies. 3 An 
3 For a precise definition of a geodesically incomplete space see e.g. |32) . 
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inextendible geodesic is a geodesic parametrized by an affine parameter r such that by using up 
only a finite amount of affine parameter the geodesic extends over infinite coordinate distances. In 
a more physical language one could reformulate the above statement by saying that it takes only 
a finite amount of affine parameter r in order to reach the infinities of the incomplete space-time. 
As we will illustrate later in this chapter, the reason why the Randall-Sundrum II setup ceases to 
be geodesically complete is simply due to a specific way of gluing two patches of AdS§. One of the 
main motivations for this work was to provide an alternative to the Randall Sundrum II model that 
has the advantage of being geodesically complete while conserving the pleasant phenomenological 
features of the latter. We divide the discussion of geodesies in two parts: in subsection 13. II we illus- 
trate the effects of incomplete geodesies in the Randall-Sundrum II setup for timelike and lightlike 
geodesies. In the following subsection 13.21 we demonstrate why our setup is geodesically complete 
by looking at corresponding geodesies. Finally, we complement the discussions by illustrating the 
physics with the use of the Penrose-diagram of (the universal covering space-time of) AdS§. 



3.1 Geodesies in the Randall-Sundrum II setup 

Our discussion of geodesies in this chapter is in no sense meant to be complete. Without going 
into the details of the computations we merely intend to present the solutions of the geodesic 
equations in certain cases. For more general and more complete discussions of this issue we refer 
to the literature, see e.g. |17| I33j and references therein. We first consider lightlike geodesies in 
the Randall-Sundrum II background metric given in appendix ^ eq. (|A.1|) . Let us suppose that a 
photon is emitted at the brane at y = in the positive y-direction at coordinate time t = then 
reflected at y = y\ at the time t = t\ to be observed by an observer on the brane at time t = t%. 
In this situation t corresponds to the proper time of an observer on the brane at rest. A simple 
calculation reveals 

t 2 = 2t x = r (W - l) . (3.1) 



k 

A brane bound observer will therefore note that it takes an infinite time for a photon to escape to 
y = oo. However, parameterizing the same geodesic by an affine parameter r reveals the lightlike 
incompleteness of the Randall-Sundrum II space-time. Let the events of emission, refection (at 
y = j/i ) and arrival on the brane again be labeled by r = 0, r = n and r = T2, respectively. Using 
the geodesic equation 

d 2 r^ dr u dr p 

+ ax ax =Q; (32) 



dr 2 up dr dr 
r 2 = 2r x = 4 f 1 - e~ kyi ) • (3.3) 



an easy computation shows: 

To = 2ti 

[ : i> constant c is a remnant of the freedom in the choice of an affine parameter. 4 From the last 
equation we see that now 

lim n = \, (3.4) 

meaning that in order to reach infinity (y = oo) in the extra dimension it takes only a finite 
amount l/(ck) of affine parameter r, the expression of incompleteness of the Randall-Sundrum II 
space-time with respect to affinely parametrized lightlike geodesies. 

4 In general, two affine parameters A and r are related by A = cr + d, since this is the most general transformation 
leaving the geodesic equation 13.21 invariant. The choice of the origin of time to coincide with the emission of the 
photon only fixes d but does not restrict c. 
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In the case of timelike geodesies the inconsistency is even more striking. As shown in e.g. |18j . a 
massive particle starting at the brane with vanishing initial velocity travels to y = oo in finite proper 
time given by r p = 7r/(2k) while for a brane-bound observer this happens in an infinite coordinate 
time. To summarize, the Randall-Sundrum II brane-world model is geodesically incomplete both 
for null and for timelike geodesies. As we will see in the next section, the geodesic incompleteness 
is a direct consequence of the use of a particular coordinate system in AdS$ , the so-called Poincare- 
coordinate system, which covers only a part of the full AdS§ space-time. We will also see that the 
problem of incomplete geodesies is absent in the setup we propose in this paper. 



3.2 Geodesies in the background (12.11) and Penrose-diagram 

We would now like to answer similar questions to the ones considered in the previous section for the 
background ()2.1|) . For example, we would like to know what time it takes for light to travel from 
the brane (at #=0) in the ^-direction to a given point in the upper hemisphere with 9 coordinate 
#i, to be reflected and to return to the brane. As in the last section, t and r will denote the 
coordinate time (proper time of a stationary observer on the brane) and the affine parameter used 
for parameterizing the geodesies, respectively. Again we choose t = (r = 0) for the moment of 
emission, t± (ti) for the reflection and t 2 (^2) for the time where the photon returns to the brane. 
Due to the enormous hierarchy of distance scales in our model, one might wonder whether photons 
(or gravitons) are able to carry information from an arbitrary point on the brane to any other point 
on the brane connected to the first one by a null-geodesic in the extra dimension. For an observer 
on the brane such a possibility would be interpreted as 4-dimensional causality-violation. However, 
as we will see in the following, none of these possibilities exist in our model. Omitting all details 
we find 



t 2 = 2ii = AR V coth 



(t) 



arctan 



sinh( 



2 



cosh [f 



0] 



so that an observer on the brane will see that the photon reaches the "north pole" 6\ 
finite time 



t\ = Rjj coth 



UJTT 



( 2 



arctan 



sinh ^ 



UJTT 

~2~ 



Ru 



7T 



(3.5) 
vr/2 at 
(3.6) 



Note that due to the warped geometry, it is Ru entering the last relation and not R, so that even 
though the physical distance to the "north pole" is of the order of R it takes a time of the order 
of Ru for photons to reach it, excluding causality violation on the brane as discussed above. If the 
same geodesic is parametrized using an affine parameter we obtain 



r 2 = 2n 



2R 

C UJ 



(urn 
tanh — 
V 2 



sinh \oj 



0] 



cosh^ 



(3.7) 



such that the amount of affine parameter needed to reach 9\ = tt/2 starting from the brane is: 

(3. 



. 1 R . /^vr\ 1 R 
At — tanh J fc 



c UJ 



C UJ 



Here again c reflects the freedom in the choice of the affine parameter. 

The results formally resemble those of the Randall-Sundrum II case. However, the important 
difference is that in our case each geodesic can trivially be extended to arbitrary values of the affine 
parameter, a simple consequence of the compactness of our space. Once the photon reaches the point 
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9 = it continues on its geodesic, approaching the brane, entering the southern "hemisphere" , etc. 
It is clear that it needs an infinite amount of affine parameter in order to travel infinite coordinate 
distances. Therefore, the null geodesies in our setup which are the analogues of the incomplete null 
geodesies in the Randall-Sundrum II setup turn out to be perfectly complete due to the compactness 
of our space. The situation for timelike geodesies is fully analog to the case of the null geodesies. 

To end this section about the geometric properties of our model we would like to discuss the 
conformal structure of our space-time and point out differences to the Randall-Sundrum II setup. 
Let us review briefly the basic properties of AdS*, space-time to the extend that we will need it in 
the following discussion. 5 AdS$ space-time can be thought of as the hyperboloid defined by 

Xl + Xf - X\ - X\ -Xl-X\ = a 2 (3.9) 

embedded in a flat space with metric 

ds 2 = -dXl - dXl + dXf + dX\ + dXl + dXj, (3.10) 

a being the so-called AdS-r&dius. The global coordinates of AdS§ are defined by 

Xq = a cosh x cos r, X§ =a cosh x sin r, 

Xi = a sinh x cos ipi , X<i =a sinh x sin </?i cos (p2 , 

X3 = a sinh x sin </?i sin (^2 cos <^3 , X4 =a sinh x sin if\ sin <f2 sin 933, (3-11) 

where the coordinates are confined by < x, —it < r < ir, < <p\ < ir, < ip2 < vr, < (^3 < 2ir 
and r = — 7r is identified with r = ir. These coordinates cover the full hyperboloid exactly once. 
Allowing t to take values on the real line without the above identification of points gives the 
universal covering space CAdS^ of AdS$. 6 In these coordinates, the line element (j3.10j) can be 
written: 

ds 2 = a 2 (- cosh 2 x dr 2 + dx 2 + sinh 2 X dftf) . (3.12) 
Another coordinate system can be defined by 

X = — \l + u 2 (a 2 + x 2 -t 2 )] , X 5 =aut, 
2u 

X 4 = aux\ i = 1,2,3 , X 4 =— \l — u 2 (a 2 — 2 2 + i 2 )] , (3.13) 

2u 

with u > 0, t £ (—00, 00) and x l £ (—00, 00). In these Poincare coordinates the line element (|3.1U|) 
takes the form 

d ^ +u ^( d ^ + d x 2 )]. (3.14) 



ds = a 



In contrast to the global coordinates, the Poincare coordinates do not cover the whole of the AdS§ 
and CAdS§ space-times (31]. Prom (j3.14|) . after changing coordinates according to dy = —adu/u 
and rescaling t and x % by the AdS-r&dius a we recover the original Randall-Sundrum II coordinate 

5 We mainly follow |53|. 

6 Whenever we used the word AdSs so far in this paper we actually meant CAdSs. For reasons of clarity, however, 
we will brake with this common practice in the rest of this section. 

7 It is obvious that the metric I12. U for < < ty/2 (—it/ 2 < 6 < 0) reduces to the above line-element 13.12H under 
the following coordinate transformation: \ = w (f T 6), T — tu)/ [i?cosh (^y-)] , together with a = R/uj. 
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system given in (jA.l|) . The restrictions on y in the Randall-Sundrum II setup further limit the 
range covered by their coordinate system to the < u < 1 domain of the Poincare coordinates. 
Coming back to the global coordinates, we introduce p by 

TT 

tanp = sinhx with < p < — , (3.15) 

so that (|3.12|) becomes 

a 2 

ds 2 = (-dr 2 + dp 2 + sin 2 p dfi§) . (3.16) 

The Penrose-diagrams of AdS§ space-time and its universal covering space-time CAdS§ are 
shown in Fig. [3 see |35l I36j . While the AdS$ space-time contains closed timelike curves (denoted 




P = ■§ p — pb p = -f 

(<P1 = t) (^i = 0) 



Figure 2: Penrose diagram of AdS^ space and its universal covering space CAdS§ 

by 7 in the figure), its universal covering space-time CAdS§ does not. The arrows in the left 
diagram indicate that the lines r = — tt and r = tt should be identified. The symbol stands for 
the timelike surface p = tt/2 (spatial infinity). It is this surface which is responsible for the absence 
of a Cauchy-surface in AdS-sp&ce. We will concentrate in the following on the CAdS$ diagram. 
First we note that each point in the diagram corresponds to a 3-sphere. The shaded region indicates 
the patch covered by the Poincare (and Randall-Sundrum) coordinates. Note that the position of 
the Randall-Sundrum brane cannot be represented in a simple way in the Penrose diagram of 
CAdS§. The reason is that the u = const, hypersurfaces of the Poincare coordinates generate 
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a slicing of flat 4-dimensional Minkowski space-times while the points in the diagram represent 
(curved) 3-spheres. From the Penrose diagram it is immediately clear that the Randall-Sundrum-II 
space-time is geodesically incomplete. The timelike curves denoted by 7 emanating from the origin 
(p = 0, t = 0) will all eventually exit the shaded region after a finite coordinate time r. These 
geodesies appear inextendible from the point of view of the Randall- Sundrum space-time. The 
problems arise due to the arbitrary cutting of a space-time along the borders of a given coordinate 
patch which covers only a part of the initial space-time. Similar conclusions can be drawn for null 
geodesies, represented by lines making angles of 45 degrees with the vertical lines in Fig. |2J The 
vertical dotted line with topology R x S 3 at coordinate pb = arctan [sinh f^)] close to p = ir/2 
in the right diagram corresponds to the location of our curved three brane, t being proportional 
to our time coordinate t (see footnote EJ) . We note that since u > 50, this line is drawn by far too 
distant from p = ir/2 as a simple expansion shows: 

7T OJ7T _ , SuJTT . . _ s 

- - 2e~— + 0{e~— ) . (3.17) 

The slice of CAdS§ to the right of pb (pb < p < tt/2) is discarded and replaced by another copy of 
the slice to the left of pb (0 < p < pb) ■ 

From the Penrose diagram we can also deduce that our space-time is geodesically complete. 
Timelike geodesies (emanating from {p = 0, r = 0)) cross the brane pb at some point or return to 
the origin depending on the initial velocities of the particles that define them. In both cases (as in 
the case of null geodesies) there is no obstacle to extend the affine parameter to larger and larger 
values. 

As we already mentioned, the absence of a Cauchy surface in CAdS'-space-time is due to the 
existence of the timelike surface at spatial infinity. By construction, our space-time excludes 
J? so that the natural question arises whether it is legitimate to revert the above argument and 
conclude the existence of a Cauchy surface in our model. Though interesting, we do not intend to 
elaborate this question any further in this paper. 

4 Gravity localization on the brane 

There are at least two equivalent ways of addressing the problem of the localization of gravity in 
brane-world scenarios. Both ways have their advantages and disadvantages. The first way is based 
on a detailed study of the Kaluza-Klein excitations of the graviton. After integrating out the extra 
dimension(s) one obtains an effective 4 dimensional Lagrangian involving the full tower of Kaluza- 
Klein gravitons. Considering the low energy scattering process of two test particles on the brane 
via exchange of Kaluza-Klein excitations allows to relate the non-relativistic scattering amplitude 
to the static potential felt by the two test particles. Since the coupling of each individual Kaluza- 
Klein particle to matter on the brane is proportional to the value of its transverse wave-function 
at the position of the brane, this approach necessitates a proper normalization of all Kaluza-Klein 
modes. An advantage of this approach is the possibility of distinguishing between the contributions 
to the potential coming from the zero mode (Newton's law) and from the higher modes (corrections 
thereof). 

The second approach is based on a direct calculation of the graviton two-point function in the 
space-time under consideration, having the obvious advantage of bypassing all technicalities related 
to the Kaluza-Klein spectrum and the wave-function normalization. However, a physical interpre- 
tation of the effects of individual Kaluza-Klein modes from the point of view of a 4-dimensional 



Pb = arctan 
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observer is, to say the least, not straightforward. Finally, due to the equivalence of the two ap- 
proaches it is clearly possible to fill this gap by reading of the Kaluza-Klein spectrum and the 
(normalized) wave-functions from the two-point function by locating its poles and determining 
corresponding residues. 

In the setup considered in this paper we choose to work in the second approach, the direct 
evaluation of the Green's functions, due to additional difficulties arising from the non-Minkowskian 
nature of the induced metric on the brane. The lack of Poincare invariance on the brane clearly 
implies the absence of this symmetry in the effective 4-dimensional Lagrangian as well as non- 
Minkowskian dispersion relations for the Kaluza-Klein modes. 

4.1 Perturbation equations and junction conditions 

We would like to study the fluctuations of the metric in the background (|2.1|) defined by: 

ds 2 = -a 2 dt 2 + R 2 vl 2 [rjij + 2#y] W + R 2 d6 2 . (4.1) 

Hij transform as a transverse, traceless second-rank tensor with respect to coordinate transforma- 
tions on the maximally symmetric space with metric rjij, a 3-sphere in our case: 

rpHij = 0, rf^iHjk = 0, (4.2) 

where V denotes the covariant derivative associated with the metric r/y on S 3 . The gauge invariant 
symmetric tensor Hij does not couple to the vector and scalar perturbations. Its perturbation 
equation in the bulk is obtained from the transverse traceless component of the perturbed Einstein 
equations in 5 dimensions: 

2 

SRmn - -Mg M N = . (4.3) 
Specifying this equation to the case of interest of a static background we find: 

where A denotes the Laplacian on S 3 . The last equation is valid in the bulk in the absence of 
sources. We can now conveniently expand Hij in the basis of the symmetric transverse traceless 

tensor harmonics on S 3 [S3 EH1 111] : 



where the satisfy 



H tJ = yzyz^(t,e)t^, (4.5) 



1=2 A 



Af^ + k?fV X) = 0, kf = 1(1 + 2) - 2, 1 = 2, 3, ... 

rf^iff^ = 0, ififp = 0, if$ = 0. (4.6) 

The sum over A is symbolic and replaces all eigenvalues needed to describe the full degeneracy of 
the subspace of solutions for a given value of /. Introducing the expansion (|4.5|) into ()4.4|) and using 
the orthogonality relation of the tensor harmonics 



^ v ik r? l f^Tj!; x) d s <p = 8 u '6^, (4.7) 
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we obtain: 

cH it^ R 2 \a 7 / -^c/7 

This equation has to be compared to the equation of a massless scalar field in the background (|2.1|) . 
After expanding the massless scalar in the corresponding scalar harmonics on S 3 we recover (j4.8|) 
with the only difference in the eigenvalue parameters due to different spectra of the Laplacian A for 
scalars and for tensors. Motivated by this last observation and in order to avoid the technicalities 
related to the tensorial nature of the graviton Hij we confine ourselves in this paper to the study 
of the Green's functions of a massless scalar field in the background ()2.1j) . 

The differential equation (|4.8|) alone does not determine Hij uniquely. We have to impose proper 
boundary conditions for H^. While imposing square integrability will constitute one boundary 
condition at 9 = ±7r/2, the behavior of H^ on the brane will be dictated by the Israel junction 
condition [10] : 

8vr / 1, 



co - cj , (4.11) 
: 1- ( Vij + 2Hij) + H'^ : = |^c ( Vij + 2H lj ) . (4.12) 



Here and g^ u denote the extrinsic curvature and the induced metric on the brane, respectively. 
T^y is the 4-dimensional stress-energy tensor on the brane 8 and the symbol : . . . : is defined in 
()2.15|) . In our coordinate system the extrinsic curvature is given by: 

K ^n- <"»> 

Hence, the non-trivial components of the Israel condition become: 

± + 2 Hij ) + H'^ : = ^ 
Separating the background from the fluctuation in (|4,12|) we find: 

: - : = 3AP C0 > (413) 
:Hlj-. = 0, (4.14) 

where we only used the continuity of Hij on the brane. While eq. (|4.13j) directly coincides with 
eq. (|2. 12|) of section ()2.3|) . eq. (|4. llf) turns out to be a linear combination of (|2.12|) and ()2.13j) . 

Eq. ()4.14j) taken alone implies the continuity of H[j on the brane. If in addition the fluctuations 
are supposed to satisfy the 9 — > —9 symmetry, this condition reduces to a Neumann condition, as 
for example in our treatment of the scalar two-point function in the Randall-Sundrum II case (see 
appendix 0) . 

In the next section we are going to find the static Green's functions of a massless scalar field in 
the Einstein static universe background. This will serve as a preparatory step for how to handle the 
more complicated case of a non-invertible Laplacian. More importantly, it will provide the necessary 
reference needed for the interpretation of the effect of the extra dimension on the potential between 
two test masses on the brane. 



T = co, and T l j = 8' jC in the notation of section il'i.Mt . 
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4.2 Gravity in the Einstein static universe 



As announced in the previous section, we will concentrate on the massless scalar field. Our aim 
is to solve the analog of Poisson equation in the Einstein static universe background. Due to its 
topology 8 x S 3 we will encounter a difficulty related to the existence of a zero eigenvalue of the 
scalar Laplacian on S 3 necessitating the introduction of a modified Green's function. 9 Note that 
in the more physical case where the full tensor structure of the graviton is maintained no such step 
is necessary since all eigenvalues of the tensorial Laplacian are strictly negative on S 3 , see (|4.6j) . 



4.2.1 Definition of a modified Green's function 

We take the metric of Einstein's static universe in the form: 

ds 2 = g^dx^dx" = -dt 2 + A 2 dn\ (4.15) 

with dft 2 being the line element of a 3-sphere given in Q2.2|) and A being its constant radius. The 
equation of a massless scalar field then becomes: 

' d lx y^g^d v u{t : x)]=j{t,x) (4.16) 



V~9 

which in the static case reduces to 

Vu(x) = Au(x) = j(x) , (4.17) 

where A is the scalar Laplacian on S 3 and we choose to write x for the collection of the three angles 
on S 3 . All functions involved are supposed to obey periodic boundary condition on S 3 so that in 
Green's identity all boundary terms vanish: 

777/^^3^ / I Z„,f=f\ -j-> „.f^t\ j3 



g (V v{x)) u{x) d 6 x = J ^f^gv(x)Vu(x)d 6 x = J ^f^gv{x) j(x) d 3 x . (4.18) 

Here and in the following bars denote complex conjugate quantities. The operator D trivially allows 
for an eigenfunction with zero eigenvalue (the constant function on S 3 ): 

Vu (x)=0, j\f^gu Q {x)uQ{x)d 3 x = l. (4.19) 

Then, since the homogeneous equation has not only the trivial solution (u = 0) which satisfies 
the periodic boundary conditions in the angular variables, the operator T> cannot be invertible. 
Therefore, there is no solution to the equation (|4.17|) for an arbitrary source. In order to still define 
a "modified" Green's function we have to restrict the possible sources to sources that satisfy the 
following solvability condition: 

^gj(x)u (x)d 3 x = 0. (4.20) 

Formally this condition can be obtained by replacing v(x) by the non-trivial solution of the homo- 
geneous equation uq(x) in (|4.18|) . We now define the modified Green's function as a solution of 

V x Q{x : f ) = S ^Ez£l _ UQ ( f ) u (x') . (4.21) 



3 For an elementary introduction to the concept of modified Green's functions see e.g |41) . 
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Replacing now v(x) by Q(x,x') in Green's identity (|4.18f) we obtain (after complex conjugation) 
the desired integral representation for the solution of (j4.1Tf) : 



u(x) = Cuq(x) + / \f—g j(x) G(x, x) d 3 x , 



where C is given by 



Note that the source 



C 



j \f^gu(x") uo(x') d 3 x" . 



. 5 3 {x-x') 
jo{x) = = u {x)u (x ) 



(4.22) 



(4.23) 



(4.24) 



trivially satisfies the solvability condition ()4.2(jj) and represents a point source at the location x = x' 
compensated by a uniform negative mass density. To fix the normalization of uq(x) we write: 



uo(x) = N $ioo(x) with / 0?$i OO (x)$ioo(x)d 3 x = 1, 



(4.25) 



where <E>ioo = — L — and so 



l 

7T 



gu (x)u (x)d 3 x = I $ 100 {x) <S> 100 (x) d 3 x = \N \' Z A 6 = 1 



|2 ,3 



|2 .3 



(4.26) 



so that 



u (x) 



1 



A 3/2 



^ioo(x). 



(4.27) 



In order to solve the differential equation (|4.21j) defining the modified Green's function Q{x, x'), we 
expand in eigenfunctions of the Laplace operator on S 3 , the so-called scalar harmonics &\i m with 
properties (see e.g. |l2l 133] ): 



A<S> xlm = (1 - \ 2 )$xim 



A = 1, 2, . . . ; I = 0, . . . , A - 1; m = -I, . . . , I. 



(4.28) 



Our ansatz therefore reads: 



oo A— 1 I 



G(X,X) = ^^2 ^ $\l m (x)cxi m (x) 



A=l 1=0 m=-l 



Introducing this in (|4.21|) we obtain 



(4.29) 



oo A— 1 I 

EE E 

A=l 1=0 m=-l 



A 2 -l 
A 2 



C\lm(x') 



§\lm{x) 



5 3 {x — x') 



Uq(x) Uq(x') . 



(4.30) 



If we now multiply by $\'i' m '(x) y/rj and integrate over 5 3 we find 

A' 2 - 1 1 — f _ — 

J^—C\>l'm>(x') = -^$\'l'rn'(x') ~ U (x') / y/rj U (x) X >l>m> (x) d 3 X, 

S3 



(4.31) 
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where we made use of the orthogonality relation of the scalar harmonics 

Vrj$\im(x) ®\'Vm>{x) d 3 x = 5\\t8l V 5 mm >. 



(4.32) 

In the case A' = 1 (and vanishing I' and m') the above equation is identically satisfied for all values 
of cioo(^)- 10 

In the case A' 7^ 1 the coefficient cxn'm'i^) follows to be 



CX'l'm' (x 1 ) 



A(l - A' 2 ) 

so that the formal solution for the modified Green's function can be written as 

oo A— 1 I 



A=2 1=0 m=-l 



$\lm{x)<&\lm{x") 

A(l - A 2 ) 



(4.33) 



(4.34) 



Due to the maximal symmetry of the 3-sphere, the Green's function Q{x,x') can only depend on 
the geodesic distance s(x,x') £ [0, ir] between the two points x and x': 



cos s = cos 921 cos (fi + sin (pi sin (p-y cos f3 , 
cos P = cos (p2 cos if 2 + sin 992 sin ip' 2 cos(</?3 



^3) 



(4.35) 



which in the case (f2 = ^2 anc l ^3 = ^3 clearly reduces to s = ipi — tp\_. Indeed, the sum over I and 
m can be performed using 11 



A-l I 



E E ^Aim(2)^ 



\lm (*) 



A sin (As) 
2^ 



sin s 



1=0 m=-l 

such that even the remaining sum over A can be done analytically: 



Q{s) = Q{x,x! 



1 



1 



8vr 2 A Air A 



(1 ) cot s 



7T 



s G [0, vr] . 



To interpret this result, we develop G(s) around s = obtaining 
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e(s) = - A 



1 3 s , Q\ 

-4^ + 8^ + l2^ + 0(s) 



(4.36) 



(4.37) 



(4.38) 



By introducing the variable r = sA and by treating Q{s) as a gravitational potential we find 



1 dQ(s) 



1 



+ 



1 



+ 0(r/A 3 ). 



(4.39) 



A ds Airr 2 12ttA 2 
We notice that for short distances, r <C A, we find the expected flat result whereas the corrections 
to Newton's law become important at distances r of the order of A in the form of a constant 
attracting force. 13 

10 This means that we have the freedom to choose cioo(S') freely. Our choice is c-km(x') = without restricting 
generality since from eq. il 1.221 we immediately conclude that for sources satisfying l|4.20^ there will never be any 
contribution to u(x) coming from cioo(a?')- 

"See e.g. gg. 

12 Note that the constant term in the expansion l|4.88|l is specific to our choice of cioo(x). 

13 There exists numerous articles treating the gravitational potential of a point source in Einstein's static universe 
(see 1111 1451 and references therein). We only would like to point out here similarities between our result (14.3711 and 
the line element of a Schwarzschild metric in an Einstein static universe background given in (44) . 
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4.3 Modified Green's function of a massless scalar field in the background space- 
time flUU) 

We are now prepared to address the main problem of this work namely the computation of the 
modified Green's function of a massless scalar field in the background space-time (|2.1|) . Since our 
main interest focuses again on the low energy properties of the two-point function, we will limit 
ourselves to the static case. The main line of reasoning is the same as in the previous section. Due 
to the fact that our space has the global topology of a 4-sphere S 4 , we again are confronted with 
a non-invertible differential operator. We would like to solve 



Vu(ipi,9) 



+ 



1 



1 







RlliPY ' R 2 a(9)j{9) 3 89 



a{9) 1 {9) 



3 du((pi,6) 
89 



(4.40) 



where the independent angular variable ranges are < ipt < it; < ipi < it; < 993 < 2n; —ir/2 < 
9 < 7r/2. Every discussion of Green's functions is based on Green's identity relating the differential 
operator under consideration to its adjoint operator. Since T> is formally self-adjoint we have 



7T 2-7T 



dipt J dipi J dips I d9^J—g (Dv)u — v{Vu) 
000 

2tt 



/[ R? 
dipi j d(p 3 -^y/rj 

n 



+ ... 



(4.41) 



where we dropped the arguments of u and v for simplicity. The dots in (|4.41|) refer to boundary 
terms in the variables ipi and since we again employ an eigenfunction expansion in scalar harmonics 
on S 3 , these boundary terms will vanish. In order to find an integral representation of the solution 
u((pi,9) of (|4.4Uj) we have to impose appropriate boundary conditions on u and v at 6 = ±7r/2. For 
the time being we assume this to be the case such that all boundary terms in ()4.41|) vanish and 
proceed with the formal solution of (|4.40f) . We will address the issue of boundary conditions in 9 
in detail in appendices 1X1 and iBl 

In the following we will collectively use x instead of (<pi,9). The homogeneous equation 
Vu{x) = does not have a unique solution under the assumption of periodic boundary condi- 
tions. In addition to the trivial solution {u = 0) we also find 



V uq{x) = , 



-g uq(x) uq(x) d A x = 1 . 



(4.42) 



Therefore, the corresponding inhomogeneous equation ()4.40[) does not have a solution unless we 
again restrict the space of allowed sources: 



-gj(x)u (x)d 4 x = 0. 



(4.43) 



As in the last section, this condition can be obtained by replacing v by the non-trivial solution of 
the homogeneous equation uq in ()4.41|) . We now define the modified Green's function by 



V x Q(x,x) 



5 4 (x-x') 



u (x)uo(x') 



(4.44) 
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From Green's identity (|4.41[) . with v(x) given by G(x, x'), we again obtain after complex conjugation 
the desired integral representation: 



with 



As before, the source 



u(x) = Cuq(x) + / \J — g j(x')Q(x\ x) d a/, 



C 



J \/—g u(x) uo(x') d A x' . 
5 4 (x — x') 



uq(x) uq(x') 



(4.45) 



(4.46) 



(4.47) 



satisfies the solvability condition Q4.43JI by construction. The normalization of the constant mode 
uq(x) is slightly more involved than before due to the nontrivial measure a(9) , y(9) 3 in the 9 inte- 
gration. By inserting 

uo(x) = N $ m ((pi)xi(B) with X i(0) = 1 (4.48) 
in the integral in H4.42|l we obtain 



iV n 



2uj 



RRfj tanh(^f) 
For the solution of eq. (j4.44|) we use the ansatz 



(4.49) 



oo A— 1 I 



(4.50) 



A=l 1=0 m=-l 



where from now on we decide to write all arguments explicitly. After inserting this in Q4.44JI we 
find 



oo A— 1 



EE E 

A=l 1=0 m=-l 



A 2 — 1 1 



5^(^-^)5(9-9') 
= j= uo((pi,d)u [(Pi,9 )■ 



1 9 1 'a(eM9) 3 - CxMd > 9,) 



00 



®Xlm(<Pi) 

(4.51) 



After multiplication by $A'J'ra' (</>i) an d integration over we obtain 



A /2 ~l ( , „ 1 1 d 

?2 777xV> C XVm'vPii ^ 9 ) H 



R(,i{ef 



i 



i 



R 2 a(9)j(9) 3 89 
5(9 - 9')^x'i'm'( l f'i) ~ uaiv'i "' ' 



09 



(4.52) 



" R%Ra(9) 1 (9Y 
We now have to distinguish the cases A' = 1 and A' 7^ 1 



y/rj u (ifi , 9) ^x'l'm' ((Pi)d 3 ipi. 
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1. (X'l'm 1 ) = (100). 

In this case, the last term on the right hand side of Q4.52JI will give a non-vanishing contribu- 



tion: 



1 1 



a 



R 2 a(9)j(9) 3 89 



a(9)j(9) 



86 



RlRa(e) 7 (ey 



M6 - 9')$ 100 (ti)- (4.53) 



|^br*ioo(^)xiWxi(^) / v^ $ ioo(^)^ioo(^)rf <Pi ■ 



By defining g^(6,9') by the relation 



^*ioo(A (1) (m') 



RfjR 



we obtain the following differential equation for g^(9,9'): 



1 



1 8 



R? a(9)-f(9) 3 89 
where we used 



MO) 



; dg^(9,9') 
89 



-Me 



(4.54) 



xmhtf), (4.55) 



2lo 



tanh (If-) 

Note that we defined Xi(9) in such a way that 



X i(9), (xi(0) = l) 



(4.56) 



2 

j a(9) 7 (9f Xl (9) X i(9)d9 



1. 



(4.57) 



2. (X'l'm') + (100). 

Due to the orthogonality of $100 and $A'i'm' 011 5' 3 5 the l as t term on the right hand side of 
(j4.52|) vanishes. We therefore have 

A ,2 -l 



1 



cy i> m '(tp'i,9, 9') + 
1 



8 



RfjRa(9H9y 



R 2 a(0h(0) 3 89 
8(9 - tf^yvM) . 



a(e) 7 (e) 



3dc \'l' m '(y'i,9,9') 

89 



Introducing g( x '\9,9') again via 



RfjR 



we see that g^'\6,9') has to satisfy 
] - X ' 2 r 9 W(9,9>) + ' ' 



8 



r 2 a(e)-i(ef 8e 



a(e) 7 (ey 



M x '\o,e') 



86 



5(9 - 9') 
a(9) 1 (9f 



(4.58) 



(4.59) 



(4.60) 
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Combining the above results for A = 1 and A / 1 we are able to write the formal solution of 
(|4~H^) as 

-. oo A— 1 I 

QiPi^eJ) = ^EEE Ww)*AJm(A (A) (M0- (4.6i) 

U a=1 1=0 m=-l 

We obtain a further simplification of this formal solution by employing the spherical symmetry 
on see eq. ()4.36j) . leaving us with a representation of the two-point function by a Fourier sum, 
which is natural for a compact space without boundaries: 

G( S , e, e>) =9(^ e, e<) = E (e, e>) . (4.62) 

Finding the formal solution (|4.62j) was straightforward apart from minor complications inherent 
to the use of a modified Green's function. To solve the differential equations (|4.55|) and (|4,60|) by 
imposing appropriate boundary conditions again is a routine task without any conceptual difficul- 
ties, though slightly technical in nature. Replacing back this solution in (|4.62|) we are left with a 
Fourier sum which at first sight looks intractable, anticipating the fact that the solutions g^(6, 6') 
(for A 7^ 1) are given by hyper geometric functions. Nevertheless it is possible to extract the desired 
asymptotic information from the sum (|4.62|) . 

In order not to disturb the transparency and fluidity of the main article we provide large parts of 
the technical calculations in four appendices. In appendix EI we report similarities and differences 
in the evaluation of the two-point functions between our case and the case of Randall-Sundrum, 
since the latter served as a guideline for handling the more difficult case under consideration. The 
solutions of eqs. (|4.55f) and ()4.6(J|) are presented in appendix [B] an d the evaluation of the Fourier 
sum (|4.62j) at distances exceeding the size of the extra dimension in appendix O Eventually, 
appendix [D] contains the evaluation of the Fourier sum (|4.62|) for distances smaller than the extra 
dimension. In this way, we can offer the reader less interested in the details of the computations 
to have the main results at hand. 

From its definition ()4.44j) we understand that the Green's function (|4,62|) can be considered 
as the response of the scalar field to the combination of a point-like source located at coordinates 
(ip^O') and a delocalized, compensating negative contribution. Since we would like to see the 
response to a point-like particle on the brane we put 8' = in eq. (|4.62|) and explicitly write the 
A = 1 term: 

with S[s,9,uj] given by (|C.1|) of appendix [0 (see also (|B.13|) of appendix[B|). The general result for 
the sum S[s,0,uj] obtained in appendixlUlis 

Sis, 0, uA = lim Sis, 6, u>] = ; ( cos s sin s 

1 1 0^0 1 1 wl- z(0) V 2 4 

- — z(0) - 5ln[l - 2(0)1 sins - -z(0)5 limMs,0,u;] , (4.64) 
where we have z(Q) = tanh 2 (^) and where we refer to (|(J.9|) for the definition of R[s,9,u]. The 
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first term in Q4.64[l is the zero mode contribution 14 and we see that it reproduces exactly the 4- 
dimensional static Green's function of Einstein's static universe given in (|4.37|) . The other two 
terms are the contributions from the higher Kaluza-Klein modes. 

We first concentrate on the case where s ~ 1 or what is equivalent r ~ Rjj. Since one can easily 
convince oneself that R{s, 0, cj] ~ [1 - z(0)]° in this regime, the contributions of the higher Kaluza- 
Klein modes are strongly suppressed with respect to the zero mode contribution. This means that 
as in the case of the Randall-Sundrum-II model it is the zero mode which dominates the behavior of 
gravity at distances much larger than the extra dimensions R. The main difference to the Randall- 
Sundrum-II case is that the zero mode of our model not only gives rise to the typical 4-dimensional 
1 /r singularity but also accounts for the compactness of space by reproducing the Einstein static 

c 



universe behavior (|4,37|) , This result is somewhat surprising given the extreme anisotropy of our 
manifold. Due to the fact that the distances between two arbitrary points on the brane are of the 
order of R, one might intuitively expect that the extra dimension can be effective in determining 
also the large distance behavior of gravity (on the brane). As we could show by direct calculation 
the above expectation turns out to be incorrect. 

Next we consider physical distances r much larger than the extra dimension r 3> R and much 
smaller than the observable universe r <C Rjj in which case the results for R{s, 0,u] can be seen to 
be: 



R[s,0,u] 



7T 7T 1 — z(0) 



2s< 
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+ o 



[1 - z(0)f 
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8 - 6 In 2 - 6 In 

S [l-z(0)]- 1/2 



s[l-z(0)] 



-1/2 



(4.65) 



1 /2 

valid for [1 — z(0)\ ' <s<l. The zero mode contribution (to S[s, 0,oj]) in this regime is simply 
the constant obtained by setting s = in the first term of (|4,64j) so that we obtain: 



S[s, 0, uj] 



7T 0(0)2 

ul-z(0) 



l + i + i[4-31n2-31ns] +0 ' " 

2s z 



(4.66) 



where we introduced s = s [1 — z(0)] . Inserting this result in (|4.63|) and using physical distance 
r = Rjjs instead of s we obtain 
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tanh («f) 



(4.67) 



14 At this point we have to explain what we mean by "zero mode" in this context, since due to the asymmetrically 
warped geometry, the spectrum of Kaluza-Klein excitations will not be Lorentz-invariant and strictly speaking, 
different excitations cannot be characterized by different 4-dimensional masses. More accurate would be to say that 
for a given value of the momentum eigenvalue A there exists a tower of corresponding Kaluza-Klein excitations with 
energies given by with n = 0, 1, . . .. In our use of language the "zero mode branch" of the spectrum or simply the 
"zero mode" is defined to be the collection of the lowest energy excitations corresponding to all possible values of A, 
that is by the set { (Eq, A) , A = 1, 2, . . .}. The definition is readily generalized to higher branches of the spectrum. 



20 



Since we would like to compare our result with the corresponding correction in the Randall- Sundrum 
II case, we introduced the dimensionless distance variable r = rui/R, the physical distance measured 
in units of the AdS-radius, in the last line of the above result. We see that in complete agreement 
with the Randall-Sundrum II scenario, our setup reproduces 4-dimensional gravity at large distances 
with extremely suppressed corrections. The only remnant effect from the different global topology 
manifest itself through the factors of tanh (^j-) and coth (^) which are very close to 1. We 
furthermore emphasize that apart from these deviations, the asymptotic we obtained coincides 
exactly with the asymptotic for the case of a massless scalar field in the Randall-Sundrum-II 
background, see (|A.20j) and e.g. gHl E3 EH1 SHI • From the factor wcoth (^f) /R in eq. (|4~F7)) we 
see that also the relation between the fundamental scale M and the Planck-scale Mpi gets modified 
only by the same factor of tanh (^y-): 

fi, = M 6 — tanh 



M l Pl = M 6 - tanh — , (4.68) 
lo V 2 / 

where we remind that R/lo is nothing but the AdS-Radius. 

Eventually, we treat the case of distances inferior to the extra dimension r <C R- After using 
the result (|D.12|) of appendix [D] in H4.63|) , a short calculation reveals 

GM0) = ^RRj ~ *J R*8* + W (469) 

a result that has to be compared to the characteristic solution of the Poisson equation in 4- 
dimensional flat space. In n-dimensional flat space one has: 



A 



1 




S(r), r=\yX\ (4.70) 



(n - 2)V S n-ir n - 2 _ 

with Vgn-i = 2tt~2 /T[Q] denoting the volume of the n—1 sphere. Specifying to n = 4, we recover 
the correct prefactor of — 1/47T 2 in (|4.69j) multiplying the 1/r 2 singularity. 

Finally, we mention that in none of the considered cases we payed any attention to the additive 
constant in the two-point function on the brane. The arbitrary constant entering the solution 

(8, 6') can always be chosen in such a way that gW (0, 0) vanishes on the brane (see appendix[Bj). 



5 Conclusions 

In this paper we considered a particular brane world model in 5 dimensions with the characteristic 
property that the spatial part of the space-time manifold (including the extra dimension) is compact 
and has the topology of a 4-sphere S 4 . Similar to the original Randall-Sundrum II model, the 3- 
brane is located at the boundary between two regions of AdS^ space-time. The coordinates of 
AdSz used by Randall and Sundrum are closely related to the so-called Poincare coordinates of 
AdS§. While the extra dimension in this set of coordinates provides a slicing of AdS§ along flat 
4-dimensional Minkowski sections (resulting in a flat Minkowskian induced metric on the brane), 
their disadvantage is that they do not cover the whole of AdS$ space-time. 

The coordinates we used in this paper are the global coordinates of AdS§ known to provide a 
global cover of the AdS§ space-time. In this case the "extra" dimension labels different sections 
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with intrinsic geometry KxS 3 , the geometry of Einsteins static universe. The induced metric on 
the 3-brane in our setup is therefore also given by R x S" 3 . 

As we illustrated with the use of the Penrose-diagram of AdS$, the incompleteness of the 
Poincare patch is at the origin of the incompleteness of the Randall Sundrum II space-time with 
respect to timelike and lightlike geodesies. Moreover we were able to demonstrate that the setup 
considered in this paper provides an alternative to the Randall-Sundrum II model which does not 
suffer from the drawback of being geodesically incomplete. The latter point was part of the main 
motivations for this work. 

The spatial part of our manifold is characterized by an extreme anisotropy with respect to one 
of the coordinates (the extra coordinate) accounting for thirty orders of magnitude between the 
size of the observable universe and present upper bounds for the size of extra dimensions. 

Another interesting property of our manifold related to the anisotropy of its spatial part is 
the fact that any two points on the brane are separated by a distance of the order of the size 
of the extra dimension R regardless of their distance measured by means of the induced metric 
on the brane. Despite the difference in the global topology, the properties of gravity localization 
turned out to be very similar to the Randall-Sundrum II model, though much more difficult to 
work out technically. We computed the static (modified) Green's function of a massless scalar field 
in our background and could show that in the intermediate distance regime R <C r <C Ru the 
4-dimensional Newton's law is valid for two test particles on the brane, with asymptotic corrections 
terms identical to the Randall-Sundrum II case up to tiny factors of tanh (^). We could also 
recover the characteristic 5-dimensional behavior of the Green's function for distances smaller than 
the extra dimension r <C R. Eventually we saw that in the regime of cosmic distances r ~ Rjj, 
somewhat counterintuitive given our highly anisotropic manifold, the Green's function is dominated 
by the behavior of the corresponding static (modified) Green's function in Einstein's static universe. 

In the simple setup considered in this paper the 3-brane is supposed to be motionless. In the 
light of recent progress in the study of 4-dimensional cosmic evolution induced by the motion of 
the brane in the bulk, it would be interesting to explore this possibility and see what kind of 
modifications of our results we would have to envisage. Finally a related, important question which 
would be interesting to address would be the question of stability of our setup. 

Acknowledgments: We wish to thank S. Dubovsky, P. Tinyakov and S. Khlebnikov for useful 
comments and discussions. E. R. is particularly grateful to E. Teufl for helpful advice in numerous 
mathematical questions. A. G. wishes to thank LPPC for the kind hospitality during most of 
this research. This work was supported by the Swiss Science Foundation. A. G. acknowledges 
"Fondazione A. Delia Riccia" for financial support. 



22 



A Parallels in the computation of the corrections to Newton's 
law between the case under consideration and the Randall- 
Sundrum II case 

The purpose of this appendix is to review briefly the calculations of the static two-point function 
of a scalar field in the Randall-Sundrum II background |471 149j and to compare each stage with the 
corresponding stage of calculations in the background considered in this paper. This serves mainly 
for underlining similarities and differences between the two calculations. Let us begin by writing 
down the metric of the Randall-Sundrum II model 

ds 2 = e-^r/^dx^dx" + dy 2 , -vrr c < y < Tir c . (A.l) 

Here y stands for the extra dimension while k and r]^ denote the inverse radius of AdS§ space and 
the (4-dimensional)-Minkowski metric with signature — h ++. Due to the orbifold ^-symmetry 
the allowed range of y is < y < 7rr c . 15 As it is well known, each 4-dimensional graviton mode 
in this background satisfies the equation of a massless scalar field. Therefore, for the study of the 
potential between two test masses on the brane we confine ourselves to solving the equation of a 
massless scalar field with an arbitrary time-independent source: 

d d 2 

Vu(x, y) = j(x, y) with V = e 2ky A x - 4k— + ^ • ( A - 2 ) 

Since the operator T> is formally self-adjoint, the corresponding Green's function will satisfy: 

^^/^ -w #\ 5 3 (x - x') 5(y - y') . 
Vg(x,x";y,y') = ^ Ll y ' . A.3 

We are now able to write down the usual integral representation of the solution of (|A.2|) : 



u(x,y) = J ^f^gj{x\y')G{x,x';y,y')d xdy', (A.4) 

where the y'-integration extends from to irr c and the x' integrations from — oo to oo. The 
absence of boundary terms in (|A.4|) is of course the result of an appropriate choice of boundary 
conditions for u(x,y) and G(x, x 1 ; y, y'). In the above coordinates of the Randall-Sundrum II case 
the orbifold boundary conditions together with the Israel condition imposed on the fluctuations 
of the metric give rise to a Neumann boundary condition at y = 0. One can easily convince 
oneself that in the limit r c — > oo the resulting Green's function is independent of the choice of the 
(homogeneous) boundary condition at y = ttv c . We therefore follow [6^ and use also a Neumann 
boundary condition at y = nr c . Another important point is that the Green's function we are 
considering is specific to the orbifold boundary condition and so describes the situation of a semi- 
infinite extra dimension. We decided to carry out the calculations in the semi- infinite case as 
opposed to where eventually the orbifold boundary conditions are dropped and the case of 
a fully infinite extra dimension is considered. 16 An immediate consequence of this will be that 



15 We allow for a finite r c only to impose boundary conditions in a proper way. Eventually we are interested in the 
limit r c — > oo. 

16 However as long as matter on the brane is considered the Green's function for the two setups differ only by a 
factor of 2. 
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the constant factors of the characteristic short distance singularities of the solutions of Laplace 
equations also will be modified by a factor of 2. Finally we note that this factor of 2 can be 
accounted for by an overall redefinition of the 5-dimensional Newton's constant leaving the two 
theories with equivalent predictions. We conclude the discussion of boundary conditions by noting 
that at infinity in x' we suppose that the u(x',y') and G(x,x';y,y') vanish sufficiently rapidly. 

The solution of eq. (|A.3|) can be found most conveniently by Fourier expansion in the x- 
coordinates and by direct solution of the resulting differential equation for the "transverse" Green's 
function. After performing the integrations over the angular coordinates in Fourier space we obtain: 



\X — X' 



(A.5) 



where g^ p \y,y') satisfies 



1 d 

e -4ky Q y 



4kyj?_Jp) 
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p 2 e 2k y g ^(y,y') 



5(y-y') 

p-Aky 



(0<y< irr c ) (A.6) 



together with Neumann boundary conditions at y = and y = 7rr c . The solution of (jA.6|) is 
straightforward with the general result 



, . P 2k(y>+y<) 

g (p) (y,y') = - k x 



(A.7) 
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h (f) K X (f e^-) - h (f e*-.) ^ 



where y> (y<) denote the greater (smaller) of the two numbers y and y' and Ii(z), ^(-z), K\(z), K2(z) 
are modified Bessel functions. Since we want to study sources on the brane we now set y' = and 
take the well-defined limit r c — > oo, as can easily be verified from the asymptotic behavior of I n 
and K n for large values of z. We are therefore able to write the static scalar two-point function as 



Q(x-x;y,0) 



e 2ky 
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sm(p\x-x'\)K 2 (fe^) 
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We would like to point out a particularity of the Fourier-integral (|A.8|) . By using the asymptotic 
expansions for the modified Green's functions K\(z) and K 2 {z) we have |5Uj : 



K2 (|eg 



e 2 



^p-f(^-i) 



P 



(A. 



showing that for y = the integral over p in (|A.8[) does not exist. The correct value of the Green's 
function on the brane is therefore obtained by imposing continuity at y = 0: 



Q{x - x"; 0, 0) = lim Q(x - x'; y, 0). 

y^O 



(A.10) 



Note that this subtlety is still present in our case of the Fourier sum representation of the 
modified Green's function (|4,62|) . as discussed in appendix O 
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The representation (|A.8|) is very suitable for obtaining the short distance behavior of the Green's 
function. By inserting the expansion (|A.9[) into (|A.8|) we can evaluate the integral which will give 
reasonable results for distances smaller than 1/k: 

g( a-^,o)~-^ |f _ g|2 + y2 . (A.11) 

We see that after replacing the exponential factor by 1 (valid for y <C 1/&) we recover (up to a 
factor of 2) the correct 5-dimensional behavior of the Green's function in flat 4-dimensional space 
(l4~T0l. 17 

For large distances the Fourier-representation (|A.8() is less suited for obtaining corrections to 
Newton's law. This is due to the fact that all but the first term in the expansion of K2/K1 in 
powers of p around p = lead to divergent contributions upon inserting in (|A.8j) . We therefore 
seek another method which will allow us to obtain the corrections to Newton's law by term-wise 
integration. 

The idea is to promote (jA.8|) to a contour-integral in the complex p-plane and to shift the 
contour in such a way that the trigonometric function is transformed into an exponential function. 
First we introduce dimensionless quantities by rescaling with k according to X = \x — x'\k, Y = yk, 
z = p\x — x'\: 



00 



U2 P 2Y r XT (^e Y ) 



=I[X,Y] 

Following [17j (see also [50]) we now use the relation 

K 2 [w] = K [w] + —Ki [w] (A.13) 

w 

to separate the zero mode contribution to the Green's function (Newton's law) from the contribu- 
tions coming from the higher Kaluza-Klein particles (corrections to Newton's law). Using (|A.13|) 
in I[X,Y] we find 

K (^e Y ) 2X 7sinz^i(fe y ^ 



I[X,Y] = I siaz -^L dz + ^ J —^frfdz . (A.14) 







=h[X,Y] =h[X,Y] 

The additional factor of 1/z in the integrand of I2 [X, Y] allows us to take the limit Y — > with 
the result: 

lim I 2 \X,Y] = vrX. (A. 15) 

Y^O 

However, we still need Y > for convergence in Ii[X, Y]. The next step is to use sinz = ^{e^} 
and to exchange the operation 9 with the integration over z: 18 

I«y]=9|/e»^^(fc|. (A.16) 



17 As discussed above the factor of 2 is a direct consequence of the boundary conditions corresponding to a semi- 
infinite extra dimension. 

18 This is justified since both the real and the imaginary part of the resulting integral converge. 



25 



The integrand in (|A. 16|) is a holomorphic function of z in the first quadrant (see |5(J| . p. 377 for 
details). We can therefore apply Cauchy's theorem to the contour depicted in Fig. El As we w iU 




Figure 3: Contour used in Cauchy's theorem to evaluate the integral (|A.16|) . 



show shortly, the contribution from the arc C2 = { z \ z = Re %ip , < <p < | } vanishes in the limit 
R — » 00 (the limit we are interested in). Symbolically we therefore have 
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0, 



(A.17) 



where . . . replace the integrand in (|A.16|) and the sense of integration is as indicated in the figure. 
This means that we can replace the integration over the positive real axis in (|A.16|) by an integration 
over the positive imaginary axis without changing the value of I\ [X, Y] . Substituting now z in favor 
of n according to n = —iz we obtain: 
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(A.18) 



In the first line in ()A.18|) we replaced modified Bessel functions with imaginary arguments by 
Hankel functions of real argument. In the second line we explicitly took the imaginary part after 
replacing the identities relating the Hankel functions and the Bessel functions J and Y of the first 
kind. The benefit from the rotation of the contour of integration is obvious at this stage. First, 
the last integral in (|A.18|) converges also for Y = due to the presence of the exponential function 
in the integrand. In this case we have 



h[X,0] 



2X 

7T 



n [Jim +[*($)]* 



dn . 



(A.19) 



Note that we used the Wronskian relation for Bessel functions to simplify the numerator, the general 
reference being once more |50j . Second, it is straightforward to obtain the large X asymptotic 
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(1 <C X) of the integral (|A. 19|) by a simple power series expansion in n around n = of the 
integrand (apart from the exponential factor), followed by term- wise integration. For the results 
see |4(jj . where the above integral (|A.19|) has been found in the wave- function approach: 



h[X, 0] ~ ^ + ^3 (4 - 31n2 - 31nX) + O 



hxX 
~X$ 



(A.20) 



The result for Q{X, 0,0) is therefore: 

k 2 



G(X,0,0) 



2irX 
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(A.21) 



However, we still have to verify that the contribution from the arc C2 in Fig. El vanishes in the 
limit R — > 00. Parametrizing z by z = Re l(fi , < tp < ir /2 we find: 



C 2 



K 



1 \x, 



(iRe ltp ) dip 



< R 



Q j_R(cos ip+i sin tp) 



K («^ e Y) 



dip . 



(A.22) 



Making now use of the asymptotic properties of the modified Bessel functions to estimate the 
second modulus in the last integrand in 1)A.22|) we can write: 
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(A.23) 



so that 
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(A.24) 



with C being a constant of the order of unity, independent of R and ip. To finish the estimate we 
observe that 



sin p + 
and obtain: 
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(A.26) 



This establishes the result that in the limit R — > 00 the arc C2 does not contribute to the integral 
in (|A.16|) and completes our discussion of the two-point function in the Randall Sundrum-II setup. 

We are now going to summarize briefly the complications arising when the general scheme 
of computations outlined above for the Randall-Sundrum II case are applied to the brane setup 
considered in this paper. First and foremost, the main difference is due to the fact that the induced 
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metric on our brane is not the flat Minkowski metric but the metric of the Einstein static universe. 
This implies the use of a modified Green's function. In the eigenfunction expansion, the discrete 
scalar-harmonics on S s will replace the continuous plane-wave eigenfunctions. We therefore expect 
the formal solution for the two-point function (the analog of (|A.5[) ) to take the form of a Fourier 
sum. 

The solution to the transverse Green's function (the analog of eq. (|A.6[l ) as presented in 
appendix El turns out to be again straightforward and completes the formal solution. 

When trying to distinguish between the zero mode contribution and the contribution coming 
from higher modes (in the language of the Kaluza-Klein approach) a relation similar to (jA.13|) 
proves useful. While the evaluation of the zero mode contribution poses no problems, the corrections 
coming from higher modes are much more involved this time. We face the following major technical 
difficulties: there is no analog of Cauchy's theorem in the discrete case of the Fourier sum. One 
solution to this problem is to use a variant of Euler-Maclaurin's sum rule in order to replace 
the Fourier sum by an analog Fourier integral and a remainder term also in the form of an integral. 
Now the procedure again is similar to the one outlined for the Randall-Sundrum case. Another 
technical problem, however, in connection with large distances corrections has its origin in the 
simple fact that by large distances in our setup we mean distances large with respect to the extra 
dimension R but also small with respect to the size of the observable universe Rjj. The detailed 
computations of the corrections in this distances regime can be found in appendix [O 

What concerns the computation of the two-point function at distances smaller than the extra 
dimension, there is no conceptual difference to the Randall-Sundrum case. The Fourier sum can 
be calculated analytically after inserting the corresponding large momentum asymptotic. We give 
the details in appendix iDl 
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B Solutions to the differential equations ( 14.551) and ( 14.601) 



The aim of this appendix is to obtain the solutions to the eqs. (|4.55f) and (|4.6U|) and hence to 
complete the computation of the modified Green's function. We will concentrate mainly on (|4,6U|) 
for two reasons: first and foremost because in the formal expansion (|4.62|) g^(9,9') multiplies a 
constant function in s and since we are mainly interested in the behavior of the Green's function 
on the brane, this constant is of no relevance for our considerations. Secondly because the con- 
struction of the Green's functions g^'(9,9') and g^(9,9') follows the standard procedure, so it is 
not necessary to go into details twice. We will merely state the result for g^(9,9'). The equation 
we would like to solve is 
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where A = 2, 3, 4 . . . and the independent variables 9 and 9' are restricted to the intervals —tt/2 < 
9 < tt/2, and —tt/2 < 9' < tt/2. The homogeneous equation associated with (|B.1|) can be reduced 
to the hyper geometric equation. Its general formal solution in the interval O<0<7r/2is given by 



W (A, 6) = aS + Vi(A, 9) + o| +) ^(A, 9) , 
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(B.3) 



Due to the symmetry of the background under 9 — > — 9 we can immediately write down the general 
solution of the homogeneous equation obtained from ()B.1|) in the interval —tt/2 < 9 < 0: 



^-\X,9) = a[ Vi(A, 



-9) + a { 2 \ 2 {K 



-9) 



TT 
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2 ~ 



< 0. 



(B.4) 



We now have to decide whether we want to restrict our Green's function to describe perturbations 
which also possess the symmetry 9 — > —9 of the background (as in the Randall- Sundrum case) 
or not. In the first case one can limit the solution of ()B.1|) to 9 > and continue the result 
symmetrically into 9 < 0. Note that imposing the symmetry 9 — > —9 on the perturbations goes 
hand in hand with imposing 9 — > —9 for the source which means that one should add a corresponding 
delta-function 5(9 + 9') to the right hand side of (jB.lj) . The Israel-condition (|4.14f) on the brane 

with 9 = then gives = 0, as in the Randall-Sundrum case. 

Since in our setup we do not have any convincing argument for imposing the 9 — ► —9 symmetry 
also on the metric fluctuations (the scalar field), the second case, which allows a breaking of this 



symmetry will be the case of our choice. The Israel-condition (|4.14|> now merely indicates that the 
perturbation (the scalar field) should have a continuous first derivative at the brane. 

In order to obtain a uniquely defined solution we have to impose two boundary conditions on 
g( x \9,9') one at 9 = tt/2 the other at 9 = —tt/2. It turns out that the requirement of square 
integrability of our solutions (with the correct weight-function) forces us to discard one of the two 
fundamental solutions at the boundaries 9 = ±7r/2, namely (f2(X,9). 
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Before going into details we introduce the following abbreviations: 



p(9) 



R 2 



S (9) 



R 2 



(B.5) 



Eq. (|B.1|) then becomes 



p(9)g^'(9,9')\ - (A 2 - l)s(9)g^(9,9') = 5(9 - 9') 



(B.6) 



where ' denotes derivatives with respect to 9 here and in the following. We start by taking 9' > 0. 
The case 9' < is fully analogous to the one considered up to minus signs. Our ansatz for (9, 9') 
is 
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(B.7) 



< 9 < 



2 ■ 



We impose the following boundary and matching conditions on g^(9,9') which will uniquely 
determine the coefficients A, B, C, D: 

1. continuity of g W (9, 9') at 9 = 0. 

2. continuity of g^'(9,9') at 9 = 0. 

3. continuity of g (x) (9,9') at 9 = 9'. 

4. jump condition of g^'(9,9') at 9 = 9'. 

Using the ansatz (jB.7|) we obtain from the above conditions: 
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[A(9>) + B(9>)}p {X) '(0) + C(9')^ X) >(0) 

[B(9')-D(9')]^ X) (9') + C(9')^ X) (9') 

[B(9>)-D(9')}ip{ X) '(9>) + C(9')^'(9') 



0., 
0., 
0. 

~pW) 



(Bi 



the solution of which is easily found to be: 
R 2 



A(9') 
C(9') 



^ , B(9<) = R2 



p[ X) (9') ^ A) (0)4 A) '(0) +4 A) (0)^ A) '(0) 



2 ^ A) (0Vi A) '(0) 
<pf\9>) 



(A) (A) „ 
<Pl ,<P2 >° 



cp[ X) (0) V [ X) '(0) 



-R 



(A) (A) „ 

n ,<p 2 >° 



D(9') = B(9') - R 2 



4V) 



(A) (A) n 



(B.9) 



where W 
used the relation 



> ^4)0 denotes the Wronskian of (/?^ and (p^ at 9 = 0. In obtaining HB.9j) we 



(A) (A) Q 
<Pl ^2 



(A) (A) „ 
</>i >^2 >0 



R 2 p(9) 



(B.10) 
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As already mentioned, the case 0' < can be treated in complete analogy to the case 9' > 0. 
Combining the two results gives the final expression for the solution g ( - x \9,9 / ) of (UHl) (for A = 



2,3,...): 



7 (A) (M') 



f A{\e'\)^\\e\) 

B(\e>\)^\\e\) + c(\e>\)^\\e\) 
I D(\e'\)v[ x) (\e\) 



-f < 6>sign((9 / ) < 0, 
< 0sign(0') < \9'\ 
\9'\ <0sign(0')< f. 



It is easily verified that the Green's function ()B.11|) has the property 



g W (0,9') = g 



(B.ll) 



(B.12) 



as expected for a self-adjoint boundary value problem. We will not need the result (|B.11|) in its 
full generality. We focus on the case where sources are located at the brane {9' = 0). In this case 
(IB. 1 ID reduces to: 



9 y '(0,0)-- (Xy— , (A-2,3,...). 

z <P\ (0) 

For the sake of completeness, we also give the solution g ( - 1 \9,9 / ) to (ET551) 



(B.13) 







R 2 <r(0h(0) 3 d8 



a(9) 7 (9) : 



M 1] {0,9') 



80 



a(8)j(9) 3 



5{9-9')-xi{0)h(O'). 



(B.14) 



The general procedure of finding the solution is fully analogous to the case of g( x) (9,9') (for A = 
2, 3, . . .), the only difference being that in regions where 9^9' (|B.14|) reduces to an inhomogeneous 
differential equation. In the interval < 9 < tt/2 the general solution is given by 



iP(9) = mo) + amo) + bmo), 



(B.15) 



with 



MO) 



R 2 



2wtanh(^j 



In 



M0) = l, MO) = 2ln 



cosh 
tanh 



7T 

~ ' 2 



7T 

^2 



coth^ 



7T 

- [ 2 



With these definitions, it is easy to verify that g^(9, 9') is given by: 



M®) -f < 0sign(0') < 
const. +M\0\)+M\0'\) + C { M\0\) O<0sign(0')< \9'\ 

M\0'\) \0'\<9sign(9')< § 



(B.16) 



(B.17) 



with 



C 



R 2 



W[Vi,^ 2 ,0] ' 

By W [-01 , -02 , 0] we again mean the Wronskian of ipi and ip2 evaluated at 9 = 0. 



(B.18) 
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C Detailed computation of the corrections to Newton's law 

In this appendix we give a detailed computation of the sum 



S[s,^]=]TAsm(As)^-^ 

A=2 Pi (0) 



(C.l) 



with 



, \ A-l 



[l-z(9)] 2 2Fl [^,^;\+l;z(9)] 



u(X - l)z(0)^ [1 - z(0)} 2 2 F X [A±I, A+3 ; A + 1; z 



(C.2) 



where we explicitly excluded 9 = since the Fourier sum (jC.l|) does not converge for this value 



of #. 19 We also used the definition z(9) = tanh 2 [u; 



9)]. 



Note that the sum S[s, 9, uj] contains the contribution to the two-point function coming from zero 
mode and higher Kaluza-Klein modes. Inspired by the Randall-Sundrum case we try to separate 
the two contributions by functional relations between contiguous Gauss Hypergeometric functions. 
Using (15.2.15) of reference [3U| for a = b — 1 = (A + 1) /2, c = A + 1 we obtain 



[l-z{9)\ 2 Fi 



A + 3 A + 3 



;X + l;z(9) 



+ 



A + l 
A-l 



2F1 



2F1 



A + l A + 3 



2 ' 2 
A + l A + l 



A + l 

In this way we are able to rewrite S[s,9,uj] in the form: 

S[s,9,u] = S±[s, 9,u] + S 2 [s,9,uj] , 

where 



;A + 1;*(0) 
;A + 1;^) 



Si[s,9,uj] ■ 



2 z(0) l-z{9) ^ A sin (As) 



^(fl)3 [1-Z(0)] 

1 z(0) l-z{9) 
~^z(9) 1 2 [l-z(0)] 2 f 2 A+l 



2 E ; 1 

A=2 

A sin (As) 



E 



z _(9y 

z(0)_ 



X 

2 2 



A+l A+3 . 

2 > 2 > 



X + l;z(9)] 



2F 1 [A±l ) A±3 ;A + 1;z( o)] 



2^1 [ 



A+l A+l. 

2 ' 2 



;A + l;z(0)] 



2 F 1 [A±i ) A±3 ;A + 1;z(0)] 



(C.3) 
(C.4) 

(C.5) 
(C.6) 



From the asymptotic formula ()D.1|) we infer that now the sum Si[s,9,uj] is convergent even for 
9 = due to the additional power of A in the denominator. We therefore obtain 



5i[s,0,u;] 



2 z(0)2 
wl - z(0) 



A sin (As 
^2 



A=2 



1 



2 z(0)2 
wl - z(0) 



7T — S 1 

cos s — - sin s 

2 4 



(C.7) 



Since we recognized in the last sum the two point function of Einstein's static universe, we attribute 
the contribution coming from Si[s,9,lu] to the zero mode in the Kaluza-Klein spectrum. We also 

19 Note that this mathematical delicacy about the divergence of the Green's function on the brane also applies to 
the Randall-Sundrum case of the Fourier integral representation of the two-point function. The obvious remedy is to 
keep 9 strictly greater than zero during the whole calculation and eventually take the limit 9 — > 0. 
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should mention that in S2[s,9,u] we still need 9 > for convergence since only then the factor 

A-l 

[z(9) / z(0)] 2 provides an exponential cutoff in A for the sum. 

The main challenge in the evaluation of the two-point function is therefore to tame the sum 
S2[s,0,lo]. Our strategy is the following: first we extend it from A = to oo by adding and 
subtracting the A = and A = 1 terms. Next, we replace the sum over A by the sum of two 
integrals employing a variant of the Euler-Maclaurin sum rule called the Abel-Plana formula 
(see. e.g. |51j . p. 289-290). Finally, we shall see that the resulting (exact) integral representation 
will allow us to extract the desired asymptotic of the two-point function on the brane in the distance 
regime of interest (R <C r <C Rjj)- 

As announced, we start by extending the range of the sum from to oo. Since the addend with 
A = vanishes, we only need to subtract the A = 1 term with the result: 



1 z(0)2 l-z(g) . 1 z (0) 1 - z{6) pr 
t>2 [s, 9, lo\ = — - 77^-., — In [1 — z{9)\ sin s }-- —K[s, 9, u>\ 



2uj z(9) 1-2(0) 
where R[s,6,uj] is defined by 

A sin (As 



R[s,9,u]=J2 



A=0 



A + l 



z(9) 
z(0) 



_, A 



u z (9)2i-z(oy 



2Fl [A±l ; A±l ;A + 1;z( g)] 

2-^1 [t^j A + l; 



(C.8) 



(C.9) 



We used the identity 

'A + l A + 3 



2^1 



A + l;z(0) 



[l-z(O)]- 1 2 F 1 



A+l A-l 



; A + l; 2 (0) 



(CIO) 



in the last step (see 15.3.3 of (SHi)- The next step consists of extending the sin function in ()C.9|) to 
an exponential and taking the imaginary part out of the sum (as in the Randall-Sundrum case). 
Now we make use of the Abel-Plana formula, first considering the partial sums: 



n 

lf(0,s,9,u;) + ±f(n,s,9,u;) + f f(X,s,9,u)dX 



f(iy,s,9,u) - f(n + iy,s,9,u) - f(-iy,s,9,u) + f(n - iy,s,9,u) 



where for the sake of clarity we introduced 

Xe tXs \z{9) 



f(X,s,9,uj) 



A • 1 



z(0) 



e 2iry _ i 



* 2F 1 [^,^;X + l;z(9)] 
2Fl [^,^;X + l;z(0)] 



(C.ll) 



dy 



(C.12) 



In the limit n — > oo we find: 



R[s,0,u)] 



/(A, s, 9, u)dX + 1 



f{iy,s,9,uj) - f(-iy,s,9,uj) 
e 2 ^ - 1 



dy 



(C.13) 



siii[a,0,w] 



=R 2 [s,e,uj] 
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C.l Leading order asymptotic of S {-Ri[s, 9, u>]} 

We focus first on Ri[s, 9,u] and observe that the function f(X,s,9,uj) is holomorphic (in A) in 
the first quadrant. 20 Therefore, in perfect analogy with the Randall- Sundrum case treated in 
appendix [XJ we can use Cauchy's theorem to replace the integration over the positive real axis 
by an integration over the positive imaginary axis. 21 After substituting A with ly we obtain : 



Ri[s,0, 



2/(1 ~ W) 

i + y 2 ' 



-ys 



z(9) 
z(0) 



2^1 



2 > 2 > 



l+iy;z(9) 



iFi 



dy. 



(C.14) 



The upper deformation of the path of integration is advantageous for at least two reasons. Firstly, 
through the appearance of the exponential in the integrand, convergence on the brane does no 
longer rely on 9 > and we can set equal to zero in (|C.14|) to obtain: 



i2i[s,0,u;] 



y(i - w) 
i + y 2 ' 



-ys _ 



2F1 


l+iy 
2 ■ 


1+ty. 
2 ' 


l + iy;z(0) 


2 Fx 


l+iy 
2 ' 


— 1+12/ 
2 


;l + iy,z(0) 



dy. 



(C.15) 



Secondly, the form (jC.14|) is perfectly suited for obtaining asymptotic expansions for large distances. 
Note that since s G [0, ir], we carefully avoided saying "for large s" since this would not correspond 
to the case of our interest. We try to compute the two-point function in a distance regime R <C 
r <C Rjj-i distances clearly far beyond the size of the fifth dimension but far below the size of the 
observable universe. Rewritten in the geodesic distance coordinate s this becomes u/ sinh \^f) -C 
s <C 1. It is this last relation which complicates the computation of the asymptotic considerably. 
We are forced to look at an asymptotic evaluation of ()C.15|) (or HC.14|) ) for intermediate s values 
such that we cannot make direct use of Laplace's method. The correct way of extracting the above 
described asymptotic is to expand the ratio of hypergeometric functions in (|C.15|) in a power series 
of the small parameter 1 — z(0). Using rel. (15.3.10) and (15.3.11) (with m = 1) of [HO] we find 
after some algebra 



+ 



2F1 


l+iy 
2 ' 


2F1 


l+iy - 


2 ' 


fl + iy\ 3 




2 J 


-( 


l + iy 




2 



±±M!+M; 1+^(0) 



pt;l + iy;z(0) 

2^(2) - 2V ' " + " J 
l + iy 



, 1 + ly , , . 
2 7 + 2^(^-^) +ln[l 



z(0)] 



In [1-jb(0)] 



2ip 



+ O 



In [1-2(0)] -V>(1) " 
|ln[l-z(0)] [l-z(0)] 2 } . 



2V(1 

-V(2)+^ 



1 + m/ 
2 

3 + iy 



In [1 - z(0)] 
l + iy 



[1 



(C.16) 



20 Everything but the ratio of the hypergeometric functions is clearly holomorphic. The dependence of Gauss's 
hypergeometric function on the parameters is also holomorphic so the only danger comes from the zeros of 



2-F1 l^-^-, ^i"^> A + 1; 2(0)] . These, however, turn out to be outside of the first quadrant. 
21 To be mathematically correct, one has to apply Cauchy's theorem to the quarter of the disk of radius R bounded 
by the positive real axis, the positive imaginary axis and the arc joining the points R and iR, see Fig.|^|in appendix 
1X1 The above statement will then be correct if in the limit R — > oo the contribution from the arc tends to zero. We 
will verify this in section ICJ. 41 of this appendix. 
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In the above relation 7 denotes Euler-Mascheroni's constant and ip(z) = T'[z]/T[z] the so-called 
Digamma-function. Note that we developed up to linear order in 1 — z(0) since we also want to 
compute the next to leading order term in the asymptotic later in this appendix. 

We are only interested in the imaginary part of R±[s, 0,uj] and it turns out that the first term 
in the expansion (|C16|) after insertion in (|C.15|) can be integrated analytically: 



oc 

{r? ) [ s ,om} = '*IJ 



ye ys ip 



l + iy 



dy 



ye- ys Q 



1> 



l + iy 



dy. (C.17) 



The last step is justified since also the real part of the first integral in ()C.17|) converges, as follows 
immediately from the asymptotic ip(z) ~ In z — ^ + O(-zz). Since 



1 + w. 



2z 
IT 



tanh 



(C.18) 



(jSOl, p.259, 6.3.12) we find 

00 

9{4 0) [ S ,0,u,]} = 1 J ye^ tank (^fjdy 





TT d 

~2ds 



-ys 



,0 



1 + e~*y 



it d 
2^~ds 



where we introduced the /3-function by 

P(x) = 

(see |52|, p.331, 3.311, 2.). 22 



e w 
1 + e- 



-dz 



it d 
2ds 

l) dy 



7T 
2? 



J er ys tanli(y) dy 



it d 
~2^fe 



7T TT 



1 

— +2 
s 



1 + e-*y 



dy 



(C.19) 



x + 1 



4 



(C.20) 



C.2 Leading order asymptotic of 3 {i^fs, ^, ^]} 

We can treat i?2 [s, o>] starting from (jC13|) in very much the same way as Ri[s,9,u>]. Noting that 
we can again put 6 equal to zero since the exponential factor e 2ny guarantees the convergence of 
the integral, Rq[s, 0,u;] becomes explicitly: 



R2[s,0,u] 



3 2ttj/ _ 1 



-y + w 
1 + y 2 



-ys _ 



{ -y-*y c vs 

1 + y 2 



2Fl 


"i±^,i±^; l + %y,z(p) 




2F1 


'l±pt,^y ] l +l y;z{Q) 


2F1 


2 y , 2 y ;l iy;z(0) 




2Fl 







dy. (C.21) 



2 By in l|(TT9t we mean the derivative of j3 with respect to its argument. 
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If we now employ the expansion (RTTHl) two times in (RT2T1) we find after some algebra: 



9{i4 0) [s,0,w]} 
and using ()C18|) we have 
3<[i4 0) [s,0,u;]} = 




-ys 



4 



1 + zy ^ +e y^ fl - iy 



dy 



y 



2 e 2 ^ - 1 L 



.^ s tanh(f) 



+ r" s lanli 



dy 



(C.22) 



(C.23) 



— 7T 



2ivy 



— tanh smh(ys)dy 



-7T 



-7T 





d_ 

ds 



(1 + e^y) 



s'mh(ys)dy 



cosh(ys) 

(1 + e n y) 2 



dy 



1 d 

2ds 



to 



1 d 



2ds 2 - s. 



iFi 



s s 
2,2 - -;3- -; 

7T 7T 



-1 



1-i 

U IT 

1 



i+- 

u t 



(1 + u) 



+ 



2 + A 



s s 
2,2 + -;3 + -; 

7T 7T 



While we substituted y by u according to u = e~ 7ry in the third line, we used (3.194, p. 313) of 52 
with z^ = 2, u = (3 = 1 and u = 2 =F s/tt (valid for s < 2tt) in the last line. In order to simplify 
the hypergeometric functions, we first use another relation between contiguous functions, namely 
eq. (15.2.17) of [SH] with a = 1, b = k, c = k + 1 and z = —1 and then the formulae (15.1.21) and 
(15.1.23) of [SO] together with the duplication formula for the T-function to obtain: 



k 



2 .Fi[2,fc;fc + l;-l] 



2 F 1 [l,k;k;-l] 



T[k] 



k-1 



A: 



2 F 1 [l,k;k + l;-l] 



2- fc Vvr 



r [§] r m 
l + (i-k)p(k). 



(k - 1) 0(k) 



(C.24) 



After making use of this in (|C23j) we finally obtain for the contribution to lowest order in 1 — 2(0): 



id_ r 

2ds 



7T 



7T 



(3(2~-)+ 1 + - 0(2 + 



7T 



2s z 7T 7T 



^{4 0) [s,0,u;]} 

Summarizing the results to lowest order in 1 — z(0) we therefore have 
3{<>M, W ]} = * 

9{ 4»> M ,„ 1} = i- t 

Expanding HC.26|) around s = we find: 

3{M 0| [. s ,0,c l ,]} = ^-^+O( s ), 

1 «fW^, +0W . 



7T 



(C.25) 



1 



7T 



0(2 



7T 



+ 1 + 



7T 



0(2 + 



7T 



(C.26) 



{4 0) [ S ,0,^]} 



2tt 2 7 



(C.27) 



with £ denoting Riemann's (^-function. As expected the contribution from i? 2 [s, 0, to] is sub- 
leading with respect to the one from r[ [s, 0, a;]. 
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C.3 Next to leading order asymptotic of S {Ri[s, 6, u>]} 

We would now like to obtain the next to leading order term in the asymptotic of 9 {R± [s, 6, uj}} , that 
is the term obtained from ()C.15|) by taking into account the linear contribution in 1 — z(0) of the 
expansion (|CJ. 161) . Since the calculation of the integral obtained in this way is rather cumbersome 
and anyway we do not need the full s-dependence, we do not intend to evaluate it fully Extracting 
the leading s-divergence at s = is sufficient for our purposes here. To start, we insert the linear 
term of (fTTTnj) in (jTTTSj) : 




l + iy 



, 1 - m - 
ayy— — 2 e 

1 + 2T 



-l + iy 



ys 



l + iy 



2V>(1) - 2ip 



2V>(2) - 2if> 
1 + ly 



3 + iy 



In [1 - jb(0)] - ^(1) - if>(2) + i> 



3 + ty 



ln[l - z(0)] 
l + iy 



In [1 - z(0)] 

x (C.28) 
[l-*(0)f 



After some algebra and by using the functional relation of the Digamma- function ip(z) (see |50| 

1 



tp{z + 1) = if>(z) + 



we find 



(C.29) 



3 



1 - z(0) 



ye-y^{...}dy 



(C.30) 



with 



l + iy 



+ 



+ 



i + y 2 

4 

x 

l-iy 



2ip(2) - 2ip 



2ip (1) - 2i\) 



\ + vy 

2 

1 + iy 



ln[l - z(0)] 



(1 + vy) 



In [1 - z(0)] 



In [1 - z(0)} - ip (1) - V (2) + 2ip 
l + iy 



l + iy 



2ip (1) - 2ip 



ln[l - z(0)] 



(C.31) 



We will now keep only the terms proportional to y 3 in 

OHM (y 2 m ED ) since only these will 
contribute to the leading 1/s 4 singular behavior. After a few lines of algebra we find 



3 i? 



1 - z(0) 



[1 - 2 7 - In 



[1-^(0)]] J 



y 3 e- ys % 



l + iy 



dy 



+ [1 - *(0)] J y 3 e~y s Q 



l + iy 



5? 



l + iy 



dy 



+terms involving lower powers of y. 



(C.32) 
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The first integral can be reduced to the integral (|C.17I) simply by replacing each power of y by a 
derivative with respect to s and by taking the derivatives out of the integral: 



y 3 e- ys % 



l + iy 



dy 




e~ ys Q 



l + iy 



dy 



1 2 o( S 
- + -P - 

S 7T \7T 



(C.33) 



We could not find an analytic expression for the second integral. However all we need is the first 
term of its small s asymptotic 23 and this can be easily obtained from the large y asymptotic of 
5ft [V> ((1 + vy) /2)]. Since asymptotically 



iff (z) ~ In z 



1 



oo in | arg z\ < it) 



(C.34) 



2z 12z 2 -z" 

we expect that this term will contribute logarithmic terms in s and we find after straightforward 
expansions 



l + iy 



!J 



1 



1 



(C.35) 



From the last formula ()C,35|) we learn that within our current approximation of keeping only highest 
(cubic) power terms in y, it is sufficient to take into account only the contribution coming from 
half. Therefore, we need to calculate the following integral 



| / ^e^'tmh(jfj]nydy = -~[l tanh (^) my dy] , (C.36) 



=J(s) 

where we used (|(L18|) once more. Note that the term proportional to In 2 can be accounted for by 
adding a contribution of the type of the first integral in (|C.32|) . One way to find the asymptotic of 
J{s) for small s is to integrate the following asymptotic expansion of tanhx 



tanhx = 1 - 2e~ 2x + 2e~ 4x - 2e~ &x 



+ ... = 1 + 2 £(-1) 



u=l 



which converges for all x > 0. Inserting this in the definition of J{s) we find: 



(C.37) 



J(s) 



e~ ys \ny 



l + 2^{-iye-* yu 



dy 



7 + lng | 2 y^ 1 ^+i 7 + ln(« + ^ 

c ^— ' S + Z/7T 



\TT / 7T *■ — " V H 



7T 



u=l 



In the sense s< 1. 



(C.38) 
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where we used the series expansion of the [5- function given e.g. in jS^j. The last two terms in the 
above formula are finite in the limit s — > and therefore will play no role in the asymptotic (since 
they are multiplied by a factor of 1 — z(0)). 

We are now ready to assemble all contributions to the small s asymptotic of 9 -f-R^fs, 0, 



1 — z(0) r , , , , n . IT d 3 

— jli [1 - 2 7 - In [1 - z(0)] + 2 In 2] --^ 



S 7T \7T/ 



[1 - 2(0)] 




.i±^ + 2 (7+ln ^(i + 1 ) + 2£L- i r +1 M><+J) 



7T 



l/=l 



+ o 



In s , 



8 - 6 In 2 - 6 In 



+ <?( 



Ins, 



(C.39) 



where we used the following series expansions for the /3-function and for the last term in (|C.38|) : 

/?(-) = --ln2 + 0( S ), 

\7T/ S 



IT 



0[- + l) = ln2 + 0(s) 



~ ( -ir + Mn(, + i) 4[in(2)]2 _ 7iii(2)+o(4 



l/=l 



(C.40) 



Note that only the term /3(s/7r) gave a contribution to the leading 1/s 4 divergence. Relations 
(RT261) and (IC~39l) constitute the main result of this appendix. They provide the desired asymptotic 
behavior of the two-point function in the regime of intermediate distances u>/ sinh (^p) Cs<Cl. 
A couple of remarks are in place. Clearly, (|C.26|) is not only an asymptotic result but is valid 
for all s G [0, 7r]. Secondly we were writing 0(lns/s 3 ) in ()C.39|) to indicate that we dropped all 
contributions coming from terms in the integrand (|C31I) lower than cubic order. Finally we point 
out that the development of the ratio (|C.16|) in power of 1 — z(0) indeed generates an asymptotic 
in the correct distance regime. This is clear from the fact that each new power in 1 — z(0) is paired 
with an additional power of y 2 in the integrand of 9 9, w]}. This by itself implies upon 

multiplication by e~ ys and integration over y an additional power of 1/s 2 . We therefore effectively 
generate an expansion in powers of \Jl — z(0)/s = 1/ [scosh (^j-)} ■ It is also easy to understand 
why this rough way of counting powers works. The reason is that the asymptotic properties of our 
integrals , u] | are determined only by the behavior of the corresponding integrands for 

large values of the integration variable y. The terms K [ifi ((1 + ly) /2)] and 9 [tp ((1 + ty) /2)] do 
not interfere with the above power counting since their behavior for large y is either logarithmic in 
y (for M) or constant in y (for 9). 

Finally we want to point out that the expansion in powers of 1/ [scosh (^y-)l is not useful for 
obtaining information about the Green's function for distances s -C a;/ cosh (^)- In this limit the 
asymptotic expansion of the hypergeometric functions in (|C.2j) for large values of A proves the most 
efficient way to recover the properties of a Green's function in a 4-dimensional space. See appendix 
ID] for the detailed form of the Green's function at distances smaller than the extra dimension. 
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C.4 Estimate of the arc contribution to Ri[s,6,u>] 

This final subsection is dedicated to the verification that the arc denoted C2 in Fig. |3] of appendix 
1X1 gives a vanishing contribution to the integral Ri[s,9,u] and thus justifying the representation 
(|C.15|) . We need to know the large A asymptotic behavior of the hyper geometric functions entering 
the definition of R\ [s, 8, u\ in ()C.12j) . The relevant formula (|D,1|) can be found in app endix iDl where 
we discuss the short distance behavior of the Green's function. After canceling all common factors 
from the ratio of the two hyper geometric functions we find: 



2 Ft [A±! ; A±i ;A + i ;z (0)] A + 1 



2 F 1 [Y> A f i ; A + 1 ^(°)] 



A-l 



z(0) 



X+l 
2 



1+0 T- 



1 

A 

Substituting A = Re^, we now estimate Ri[s,9,lu]: 



-v{6) 



A+1 
2 



-HO)] 



[1 



,-"(0)1 



[1 + 



-1/(0)13 



[1 



-"(0)1 2 



(C.41) 
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(C.42) 



where C is a constant independent of R and the estimate of the ratio of the two hyper geometric 
functions is valid for large enough R. We also introduced a as in appendix [D] by (|D.7|) . Note that 
a < 1 for < 8. In the last line we estimated the inverse of the square root by an arbitrary constant 
(e.g. 2), certainly good for large enough R and for all p. All what is left is to observe that 



< e 

to obtain the desired result 



In ex In Ol 
min[l, ]< sin p H cos ip for s > 0, a -1 > 1 



(C.43) 
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for R — > 00. 



(C.44) 
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D The two point function at very short distances (r <C R) 



The aim of this appendix is to calculate the behavior of the two-point function at distances smaller 
than the extra dimension r <C R- We will do so by using an appropriate asymptotic expansion for 
large parameter values of the hyper geometric functions in QC.2|) . Such an asymptotic can be found 
for example in or |54j . The general formula is 24 



2-F1 [a + n, a — c + 1 + n; a — 6 + 1 + 2n; z] 



with v defined via the relation 



2 a + b T[a-b + l + 2n] (2 



7T \ 2 



T[a- c+l + n]T[c-b + n 
{l + e^)^ (l - 



—u(n+a) 
—, X 

z (n+a) 

1 + 



(D.l) 



2(1 



Z) 2 



We need to calculate 



S[s, 9, lo] = ^ s i n (^ s ) 



A=2 



(D.2) 



(D.3) 



with the ratio ip^ (6) /<p \ (0) given in (|C.2j) . Using the above expansion (jD,l|) two times we find 
after numerous cancellations 



S[s,0,u] 
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CO 



l-z(9) \ 2 [l + e 
1 - 2(0) 



[1 



z(0)2 e 2 
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where we used the following abbreviations 

z(0) = tanh 2 



LO 



7T 



tanh 



2 

LO /IT 



2 \2 



(D.4) 

(D.5) 
(D.6) 



Note that (|D.6|) is obtained after some lines of algebra by inserting (jD.5|) in ()D.2|) . Still the sum 
in l)D.4|) diverges for 9 = but it turns out that it can be performed analytically for 9 > so that 
after the summation the limit 9^0 exists. For the sake of a lighter notation we introduce the 
symbol a by 



a 



tanh [f (j-fl)] 
tanh (if) 



(D.7) 



24 Since we noticed that the formulae given by Luke, p. 452 eq. (20)-(23) |54| and the formula given by Erdelyi 
53 differ by a factor of 2 a+b we performed a numerical check. Its outcome clearly favored Erdelyi's formula which 
we reproduced in ID. II 
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It is now straightforward to evaluate the sum over A: 



A=2 



A sin (As) 
A- 1 
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The full result can therefore be written as 
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Although the last result is given in closed form, we have to emphasize that its validity is restricted 
to distances smaller than the size of the extra dimension R. 
It is now save to take the limit 9 — > in (ID.9I). The result is 



lim Sis, 9, u] 
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4sin 2 (|) 
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(D.10) 



We note that it is the first term in the curly brackets on the right hand side of (|D.9|) which is 
responsible for reproducing the characteristic short distance singularity of the two-point function. 
The main result of this appendix is therefore: 



lim Sis, 9, u>] 

8^0 



sinh (if) 



- + ?- + 0(slns] 
s 2 



(D.ll) 



The 5-dimensional short distance singularity can even be obtained including the extra dimension. 
Apart from the logarithmic term in the curly brackets of (jD.9|) . all but the first one are finite 
at short distances (in 9 and s). Developing the denominator of this term, we recover the usual 
euclidean metric in M 4 : 
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cos s ~ - 



sinh 2 (tf) 



f + s 2 



+ . 



(D.12) 



where the dots denote higher terms in 9 and s. 
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